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Where are the roots of the Bethe Ansatz equations? 
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Changing the variables in the Bethe Ansatz Equations (bae) for the xxz six-vertex model we 
had obtained a coupled system of polynomial equations. This provided a direct link between the 
BAE deduced from the Algebraic Bethe Ansatz (aba) and the BAE arising from the Coordinate 
Bethe Ansatz (cba). For two magnon states this polynomial system could be decoupled and the 
solutions given in terms of the roots of some self-inversive polynomials. From theorems concerning 
the distribution of the roots of self-inversive polynomials we made a thorough analysis of the two 
magnon states, which allowed us to find the location and multiplicity of the Bethe roots in the 
complex plane, to discuss the completeness and singularities of Bethe’s equations, the ill-founded 
string-hypothesis concerning the location of their roots, as well as to find an interesting connection 
between the BAE with Salem’s polynomials. 

PACS numbers: 02.10.De, 02.30.1k, 05.50.+q 


Exactly integrable models provide benchmarks for dif¬ 
ferent areas of physics as statistical mechanics [ij, con¬ 
densed matteythysics [H, quantum field theory [ 3 , nu¬ 
clear physics Q, atomic-molecular physics Q and more 
recently for high energy physics through the gauge the¬ 
ory, string theory and super-Yang-Mills theories [Q. An 
important tool is the algebraic Bethe Ansatz Q culmi¬ 
nating in the Bethe Ansatz equations ii- Yet analytic 
results have been unable to scale the unsurmountable 
wall of find their roots: these have treated mostly by 
numerical methods, which are in general hard to imple¬ 
ment. 

Here we analyze the solutions of the Bethe equations 
using some theorems regarding self-inversive polynomials 
in order to answer the question made in the title. 

The BA equations for the xxz six-vertex model on a 
L X L lattice, as deduced from the ABA are 

/ sinh(Ai ^ _ yr sinh(Ai - Xk + rj) 

V sinh Xi ) sinh(Ai - Xk - 0 ) ’ 

i=l,2,...,N. (1) 

The solutions {Ai, A 2 ,..., Aat} of (U) will furnish all 2^ 
states of the transfer matrix for a lattice of L columns. 

Multiplying side-by-side these equations we get 


1^ / sinh(A^ -b ??) A ^ ^ 
M V sinhA* ) 


( 2 ) 


which suggests the following changing of variables. 


Now, Eq. ([3]) can be easily solved for the rapidities Ai, 


, , / smh -q \ 

Xi = arctanh -^— 

\Ci — cosh 77 / 



( 5 ) 


where the functions arctanh(z) and ln(z) must be re¬ 
garded as multivalued complex functions, each branch 
differing by multiples of in. The new variables Ci should 
yet to be determined. To this end we insert (O back on 
o, which gives us a system of N polynomial equations. 


L _ /_i xAT-l TT CjCfc — 2Aci -I- 1 
* IV ■ ~ 2^Cfe -I-1 ’ 

A = cosh? 7 , i = l,2,..,N. 


( 6 ) 


Here we observe that writing Ci = exp(A:i), where ki 
are the Bethe’s momenta, leaves dH]) exactly equal to the 
BAE derived in the CBA for the xxz six-vertex model [^. 
Therefore, the relations establish a direct link be¬ 
tween the Bethe states of the aba and the Bethe wave- 
functions of the CBA. Thus, all the results about com¬ 
pleteness, singularities etc. which are valid for CBA, as 
obtained by Baxter in , will be also valid for the al¬ 
gebraic version. 

For N = 1 the equations dH]) reduce to cf = 1 and ci is 
one of the L roots of unity. This means that in a periodic 
lattice of L columns, the free pseudo-particle (magnon) 
has L different rapidities given by the Bethe roots ©■ 
For N = 2 we have three coupled equations for ci and 
C2, 


sinh(Ai -b 77 ) 
sinh Ai 


i = 1,2, 




( 3 ) 


so that the Ci should be subject to the constraint equation 


Cl (ciC 2 — 2 Ac 2 + 1) + C 1 C 2 — 2Aci -1-1 = 0, 

C 2 (ciC 2 - 2Aci + 1) + C 1 C 2 - 2 Ac 2 -1-1 = 0, 

cic2 = 1 - ( 7 ) 


cf C 2 • • • — 1, (4) 

which reflects the translational invariance of the periodic 
lattice. 


From the constraint equation we can set C 2 = Wa/ci and 
the system becomes reduced to 

Pa{ci) = 0, Pa(ci)/cf = 0, ( 8 ) 
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where Pa(ci) are the following L-degree polynomials, 

Pa{ci) = (1 + UJa)Ci — 2ALOaCi ^ — 2Aci + (1 + Wa), 
Wa = e^“, a = l,2,...,L. (9) 

Notice that Pa(ci) satisfies Pa(ci) = cfPa(wa/ci), the 
pair (ci,C2) and (c2,ci) therefore representing the same 
solution of 0 . Moreover Po(ci) is a self-inversive poly¬ 
nomial since they satisfy Pa{ci) = uiacf Pa{^/ci), where 
the bar means complex-conjugation. 

It seems that explicit solutions (in terms of radicals) of 
the equation Paifii) = 0 can be written only for small val¬ 
ues of L, or for very special values of a and A. The roots 
of Pa{ci) of course can be easily obtained numerically. 
It is the self-inversive property of the polynomials Pa{ci) 
that provides crucial informations about the location of 
the Bethe roots. 

We know from the theory of self-inversive polynomials 
that their roots are all symmetric with respect to 
the complex unit circle U = {z & C : \z\ = \}. For the 
polynomials Pa(ci), however, we get a better scenario, 
since their roots can be distributed only into two ways, 
namely, either they are all on U or only two dual roots 
s and LOa/ s are not on U. In fact, the exact behavior of 
the roots depends on the values of A from which two are 
critical, 


Ai 


1 + LOa 
2 


and A 2 


L 1 LOa 
L-2 2 


( 10 ) 


and this behavior can be summarized as follows, 

i) If |A| < Ai then all roots of Po(ci) are on U and 
they are simple; 

ii) If |A| > A 2 then only two dual roots of Po(ci) are 
out U and all roots are simple; 

in) If |A| = Ai then all roots of Pa(ci) are on U but 
two roots on U can be coincident; 

iv) If |A| = A 2 then two roots of Pa(ci) are out U but 
three roots on U can be coincident; 

v) If Ai < |A| < A 2 The roots are distributed as in i) 
or ii), depending on the values of L and a. 

The proof of the statements i) and i£) can be obtained 
from theorems presented by Vieira in [l2j| , which are gen¬ 
eralizations of theorems presented previously by Lakatos 
and Losonczi in [l^ : Hi) and iv) can be verified directly 
from the factorization of Pq(ci) and v) was verified nu¬ 
merically only (since the theorems mentioned above does 
not hold in this case). The exact behavior in the case v) 
is the following: all roots of Pa(ci) will be in U if, for even 
L, a is odd and for odd L, if a is odd and a < (L — l)/2 
or when a is even and a > {L + l)/2; for a = {L — l)/2 
or a = {L + l)/2 the roots are all in P if L = 3 (mod 4); 
otherwise the roots are in U, except for two dual roots s 
and LOa/s. 

Notice that in the case in the case i), where the roots 
of Pa(ci) are all in U, if the coefficients of Pa(ci) are 
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FIG. 1. Roots of Pl{ci) and Pi(ci) for L = 20 and A = 2. 
The first figure shows that Pl{c\) is a Salem’s polynomial. 
In the second figure we have Pi(ci) behaves as a “rotated” 
Salem’s polynomial. 


integers then follows from a celebrated theorem of Kro- 
necker 1^ that all roots of Pa{ci) are indeed roots of unit 
and therefore they can be expressed by radicals. More¬ 
over, in ii) the polynomials Pa(ci) have only two dual 
roots out of U. In algebraic number theory a polynomial 
(with integer coefficients) whose roots are all on the com¬ 
plex unit circle except for two positive reciprocal roots r 
and 1/r is named a Salem polynomial. Therefore the 
polynomials Pl (ci) are Salem’s polynomials when A is a 
positive integer greater than 1. See Fig. [Tl[2]and[3l This 
is an interesting and non expected relation between the 
BAE for the two magnon state and Salem’s polynomials, 
since they are found only in a few fields of mathemat¬ 
ical physics, for instance, in Coxeter systems and the 
(—2,3, 7)-pretzel knot theory 151. 


The analysis above gives us the distribution of the 
roots of Pa(ci). The correspondent location of the Bethe 
roots Ai is provided by the formula ([5]), which can be 
seen as a conformal mapping from the complex variables 
Ci to Ai. In fact, if |A| > 1, so that 77 is real, this map¬ 
ping will send the complex unit circle U into the vertical 
line X = —I 77 I/ 2 . On the other hand, if |A| < 1 then 
([5|) will map U into a horizontal lines y = —I 77 I /2 and 
y = —|77 + *7r|/2. Notice that usually the Bethe roots are 
expected to be arranged (in the thermodynamical limit 
L ^ 1) into groups of the same real part called strings 
[ 1 ^ . A n-length string is a group of n Bethe roots all of 
them with the same real part. From what was said above 
we can see that, for |A| < 1, the roots of Pa(ci) which lie 
on U will lead to 1-string solutions, while the two dual 
roots outside U will lead to 2-string solutions. The same 
will be true for the case |A| > 1 provided we redefine a 
n-length string as group of n Bethe roots all of them with 
the same imaginary part. See Fig. |4] and [H Notice as 
well that several violations of the string hypothesis were 
already reported [l6| . This not happens in our approach, 
since the string hypothesis is not used in our deduction. 


Here we remark the importance of this analysis for the 
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• a = 1 % a = 2 #3 = 3 #3 = 4 •3 = 5 •a = 6 

FIG. 2. The roots of Pa(ci) for L = 6 and A = 2/3. For 
a = {1, 5, 6} we have |A| < Ai, so the roots are all in U. For 
a = {2,4} we have Ai < |A| < A 2 and since a is even we 
have in each case two dual roots out of U. For a = 3 we have 
I A| > A 2 , then two roots are also out of U (these roots are 0, 
and — 00 ). 



• a = 1 *3 = 2 *3 = 3 *3 = 4 *3 = 5 •a = 6 


FIG. 3. The roots of Pa{ci) for L = 6 and A = 3/2. For all 
values of a we have |A| > A 2 (the equality occurs only for 
P^ici), which does not have three roots coincident). Therefore 
all Pa{ci) have two dual roots out of U. 


range of A in the study of the completeness of the Bethe 
states and in the string hypothesis. 

Let us first to consider |A| ^ Ai and |A| 7^ A2. 
For odd L we have that Pa(ci) factors to Pa(ci) = 
(ci + uJb)Qaici), where Qaici) are L self-inversive poly¬ 
nomials of L — 1 degree and 2 b = a (mod L). However, 
the solutions ci = —wt, lead us to Ai = A2, that is, to L 
states not belonging to the two magnon sector. Hence, 
the wanted solutions are the L — 1 roots of each poly¬ 
nomial Qa{ci). A half of these solutions are related by 
permutations, so we get the exact number of solutions for 
the two magnon sector for odd L, namely, L{L — l)/ 2 . 

This same number of physical states for L even (L = 


2 k) is obtained by a more elaborate sum of terms. This 
happens because the case Wa = — 1 is special, since Pa(ci) 
collapses in this case to 

Pi/2(ci) = 2Aci(cf-2-l). (11) 


The possibility Ci = 0 leads to C2 = — 00 and the system 
dS]) is not actually satisfied. However, this would lead 
to the finite Bethe roots (Ai,A2) = {—r], 0 ) which is in 
fact a solution of the Bethe equations. The explanation 
is that we had implicitly assumed Cj 7^ 0, 1 < j < L, 
in deriving Q and thus these cases must be analyzed 
separately. Thanks to the conjugation property, we have 
as well the solution (Ai,A2) = (0,-77), so we get two 
additional solutions of the bae which must be taken into 
account. 

For the other case, Ci must be a (L — 2 )-root of unity 
and we get 


Ai = -ln 


0 - 

0 - e" 


A2 = -ln 




1 


+ 1 


( 12 ) 


where Q = for / = 1,2, ...L — 2. In this case we 

have A 2 7 ^ Ai except when j = (fc—1)/2 or j = 3(A: —1)/2, 
which can only happen if k is odd. Therefore, the case 
a = k furnishes 2(k — 1 ) solutions if k is even, but only 
2{k — 2) solutions if k is odd. Now let us consider the 
cases where a ^ k. Here follows that the equations ([9]) 
factor to Pa{ci) = — uJa)Qaici) when a is odd, but 

it not factor when k is even. For even k we have k — 1 
cases where a is even and k cases where a is odd, which 
gives us {k — 1) {2k) -I- (fc) {2k — 2) solutions. For odd k 
we have k cases where a is even and k — 1 cases where a 
is odd, which gives us {k) {2k) + {k — l){2k — 2) solutions. 
Taking into account the solutions from a = k and that 
half of the solutions are related by permutations we get, 
in both cases, a total of k{2k — 1) = L{L— l)/2 solutions, 
as expected. 

Here we remark that for L = 2, Pi{ci) is identically 
satisfied and we have a physical state with a free param¬ 
eter (see [ 1 ^ for more details). Moreover, it turns out 
that P 2 {ci) has only the singular solutions ci = as 
the solutions studied in 0 - 

Now let us consider the cases where multiple roots ap¬ 
pears. If |A| = Ai, here we have that, for L = 2k — 1, 
{k > 1), there are k — 1 polynomials (a = 2,4, ...,fc — 
2, fc -|- 1, fc -b 3,..., L) for k even and k polynomials {a = 
2,4, ..., k — 1, k, k + 2, k + 4,..., L) for k odd which factor 
in the form Pa(ci) = (1 + u}a){ci -b a;b)(ci - u}b)^Qa{ci), 
with 2b = a (mod L). For L = 2k, {k > 1), Pa{ci) 
factors as before but without double roots if a is odd. 
It turns out that for a even it now factors to Pa{ci) = 
(1 -b uja){ci — uif)'^Ra{ci) and we have k — 1 polynomials 
of this type if k is even and k polynomials if k is odd, all 
with two coincident roots identified by the indices a = 2p 
and f = p for a < k and f = k + p for a > k. Summariz¬ 
ing, the number of states for two magnon in a lattice of 
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FIG. 4. Bethe’s Roots for L = Q and A = 2/3. Notice 
that the roots of Pa{ci) which lie on U are mapped on the 
lines y = —\r]\/2 and y = —\y + i-K\/2. This leads to 1-string 
solutions. The dual roots of Pa{ci) which are out of U are 
led to pairs with the same real part, so they are related to 
2-string solutions. 


L columns is reduced to L{L — l)/2 — a where a = k — \ 
for k even and a = k for k odd when L = 2k — 1 and 
a = 2k — 1 when L = 2k. A similar count holds when 
|A| = A 2 . In this case however the number of states is 
reduced to L{L — l)/2 — /3 where /? = fc for fc even and 
/? = /c — 1 for k odd when L = 2k — 1 and f3 = 2k for k 
even and /3 = 2{k — 1) for k odd when L = 2k. 

The presence of these multiple roots in the two magnon 
sector mask the completeness [l^. However, we can use 
the 5'-matrix language in order to understand these new 
states not as two magnon states, but as free bound states 
of two pseudo-particles with the same or parallel rapidi¬ 
ties (cf = C 2 = 1). This is another physical problem. 

Finally, a short comment about the case N = 3. Now 
we have three coupled equations in ci, C 2 and C 3 plus 
the constraint equation cfc 2 C 3 = 1. In particular, self- 
inversive polynomials with Wa = e = ±I can be obtained 
by setting C 3 = 1 and C 2 = e/ci and it has the form 

P£,(ci) = - (3eA -l)cf + eA {2eA - 1 ) 

- A(2eA-l)c?-b (3eA-l)ci-£ = 0. (13) 

In fact, we have verified that all roots of Pl{ci) are in 
t/ if 0 < A < 1, but for other intervals we found that 
Pi(ci) may have two or four roots out of U. Moreover, 
by removing multiple roots of these polynomials (for each 
value of e) we got {L — 1)(L — 2)/6 states for L > 6 
and taking into account the other values of Wq the total 
is L{L — 1){L — 2)/6 states. We will present a more 
detailed study of the N >3 case in a forthcoming study. 
This study can be also generalized to more general vertex 
models as the eight-vertex model and those requiring the 
nested bae. 

{Acknowledgments). It is ALS’s pleasure to thank pro¬ 
fessor Roland Koberle for their help and advice in prepar¬ 
ing this article. The work of RSV has been supported 
by Sao Paulo Research Foundation (FAPESP), grant 
:?(b2012/02144-7. ALS also thanks Brazilian Research 



• a=1 • a=2 • a=3 


a = 5 • a = 6 


FIG. 5. Betlie’s Roots for L — 6 and A = 3/2. Now the roots 
of Pa{ci) which lie on U are mapped on the line x = —\r]\/2. 
This leads to l-string solutions. The dual roots of Pa{ci) 
which are out of U are led to pairs with the same imaginary 
parts. This leads to 2-string solutions. 
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